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Data Processing Method, Equalizer and Receiver 
Field 

[0001] The invention relates to a data processing method in an equalizer of a 
receiver, an equalizer and a receiver. 

Background 

[0002] Transmissions through multi-path channels usually experience inter- 
symbol interference at a receiver. One of the methods to mitigate the interfer- 
ence is to utilise an equalizer at the receiver to compensate for channel distor- 
tion. 

[0003] Iterative turbo methods have been used for channel equalization where a 
channel equalizer and a channel decoder exchange information iteratively. One 
of the methods consists of a soft interference canceller followed by a Minimum 
mean-square error MMSE filter (SC/MMSE) optimised with the channel de- 
coder feedback, channel response and noise level. 

[0004] The dominating computational complexity of the SC/MMSE algorithm 
lies in the computation of the inverse of the interference covariance matrix. The 
inverse is calculated for each transmitted symbol in the received signal. 

Brief description of the invention 

[0005] According to an aspect of the invention, there is provided a data proc- 
essing method in a channel equalizer of a receiver, comprising: estimating inter- 
ference from a received signal at a first observation time, creating a first covari- 
ance matrix on the basis of the estimation and defining an inverse matrix of the 
first covariance matrix and a Cholesky decomposition matrix thereof, removing 
selected covariance components from the Cholesky decomposition matrix, com- 
puting the inverse of a sub-matrix, which represents the common part of the first 
covariance matrix, and a second covariance matrix, which includes covariance 
estimates of a second observation time, by using the aid of the Cholesky 
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decomposition of the inverse matrix of the first covariance matrix, estimating 
interference from a received signal at a second observation time and determin- 
ing additional covariance components on the basis of the estimation, creating 
the Cholesky decomposition of the inverse matrix of the second covariance ma- 
trix by using unitary rotations, generating an output value of the channel equal- 
izer by utilizing information obtained with the aid of the Cholesky decomposi- 
tion of the inverse matrix of the second covariance matrix. 

[0006] According to another aspect of the invention, there is provided an equal- 
izer comprising: means for estimating interference from a received signal at a 
first observation time, creating a first covariance matrix on the basis of the esti- 
mation and defining an inverse matrix of the first covariance matrix and a Cho- 
lesky decomposition matrix thereof, means for removing selected covariance 
components from the Cholesky decomposition matrix, means for computing the 
inverse of a sub-matrix, which represents the common part of the first covari- 
ance matrix, and a second covariance matrix, which includes covariance esti- 
mates of a second observation time, by using the aid of the Cholesky decompo- 
sition of the inverse matrix of the first covariance matrix, means for estimating 
interference from a received signal at a second observation time and determin- 
ing additional covariance components on the basis of the estimation, means for 
creating the Cholesky decomposition of the inverse matrix of the second covari- 
ance matrix by using unitary rotations, means for generating an output value of 
the channel equalizer by utilizing information obtained with the aid of the Cho- 
lesky decomposition of the inverse matrix of the second covariance matrix. 
[0007] According to another aspect of the invention, there is provided a re- 
ceiver comprising: means for estimating interference from a received signal at a 
first observation time, creating a first covariance matrix on the basis of the esti- 
mation and defining an inverse matrix of the first covariance matrix and a Cho- 
lesky decomposition matrix thereof, means for removing selected covariance 
components from the Cholesky decomposition matrix, means for computing the 
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inverse of a sub-matrix, which represents the common part of the first covari- 
ance matrix, and a second covariance matrix, which includes covariance esti- 
mates of a second observation time, by using the aid of the Cholesky decompo- 
sition of the inverse matrix of the first covariance matrix, means for estimating 
interference from a received signal at a second observation time and determin- 
ing additional covariance components on the basis of the estimation, means for 
creating the Cholesky decomposition of the inverse matrix of the second covari- 
ance matrix by using unitary rotations, means for generating an output value of 
the channel equalizer by utilizing information obtained with the aid of the Cho- 
lesky decomposition of the inverse matrix of the second covariance matrix. 
[0008] According to another aspect of the invention, there is provided an equal- 
izer configured to: estimate interference from a received signal at a first obser- 
vation time, creating a first covariance matrix on the basis of the estimation and 
defining an inverse matrix of the first covariance matrix and a Cholesky decom- 
position matrix thereof, remove selected covariance components from the Cho- 
lesky decomposition matrix, compute the inverse of a sub-matrix, which repre- 
sents the common part of the first covariance matrix, and a second covariance 
matrix, which includes covariance estimates of a second observation time, by 
using the aid of the Cholesky decomposition of the inverse matrix of the first 
covariance matrix, estimate interference from a received signal at a second ob- 
servation time and determining additional covariance components on the basis 
of the estimation, create the Cholesky decomposition of the inverse matrix of 
the second covariance matrix by using unitary rotations, generate an output 
value of the channel equalizer by utilizing information obtained with the aid of 
the Cholesky decomposition of the inverse matrix of the second covariance ma- 
trix. 

[0009] According to another aspect of the invention, there is provided a re- 
ceiver configured to: estimate interference from a received signal at a first ob- 
servation time, creating a first covariance matrix on the basis of the estimation 
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and defining an inverse matrix of the first covariance matrix and a Cholesky 
decomposition matrix thereof, remove selected covariance components from the 
Cholesky decomposition matrix, compute the inverse of a sub-matrix, which 
represents the common part of the first covariance matrix, and a second covari- 
ance matrix, which includes covariance estimates of a second observation time, 
by using the aid of the Cholesky decomposition of the inverse matrix of the first 
covariance matrix, estimate interference from a received signal at a second ob- 
servation time and determining additional covariance components on the basis 
of the estimation, create the Cholesky decomposition of the inverse matrix of 
the second covariance matrix by using unitary rotations, generate an output 
value of the channel equalizer by utilizing information obtained with the aid of 
the Cholesky decomposition of the inverse matrix of the second covariance ma- 
trix. 

[0010] Embodiments of the invention are described in the dependent claims. 
[0011] The method and system of the invention provide several advantages. An 
embodiment of the invention offers an efficient method for computing a matrix 
inversion in an equalizer. The embodiment can be implemented by a pipelined, 
parallel arrangement with, for instance, one or more systolic arrays. 

List of drawings 

[0012] In the following, the invention will be described in greater detail with 
reference to the embodiments and the accompanying drawings, in which 
[0013] Figure 1 shows an example of a communication system, 
[0014] Figure 2 is a flow chart, 

[0015] Figure 3 shows an example of the SC/MMSE turbo receiver structure, 
[0016] Figure 4 illustrates a sub-matrix including the common part of two con- 
secutive covariance matrices, 

[0017] Figure 5 is a block diagram of an equalizer and 
[0018] Figure 6 is a block diagram of a receiver. 
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Description of embodiments 

[0019] With reference to Figure 1, we examine an example of a data transmis- 
sion system in which the embodiments of the invention can be applied. The pre- 
sent invention can be applied in various wireless communication systems. One 
example of such communication system is a UMTS (Universal Mobile Tele- 
communications System) radio access network. It is a radio access network, 
which includes WCDMA (wideband code division multiple access) technology 
and can also offer real-time IP-based (Internet Protocol) services, which are, for 
example, IP telephony (IPT), IP-based streaming and IP-multimedia via a con- 
nection to an IP-network. The embodiments are not, however, restricted to the 
system given as an example but a person skilled in the art may apply the solu- 
tion in other radio systems provided with the necessary properties. 
[0020] It is clear to a person skilled in the art that the method according to the 
invention can be applied to systems utilizing different modulation methods or 
air interface standards. 

[0021] Figure 1 is a simplified illustration of a digital data transmission system 
to which the solution according to the invention is applicable. This is a part of a 
cellular radio system, which comprises a base station (or a node B) 100, which 
has bi-directional radio links 102 and 104 to subscriber terminals 106 and 108. 
The subscriber terminals may be fixed, vehicle-mounted or portable. The base 
station includes transmitters, for instance. From the transceivers of the base sta- 
tion there is a connection to an antenna unit, which establishes the bi-directional 
radio links to the subscriber terminal. The base station is further connected to a 
radio network controller 110, which transmits the connections of the terminals 
to the other parts of the network. The base station controller controls in a cen- 
tralized manner several base stations connected to it. 

[0022] The cellular radio system can also communicate with other networks, 
such as a public switched telephone network or the Internet. 
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[0023] It is possible to apply the turbo principle to many detection/decoding 
problems such as channel estimation, channel equalization, detection of coded 
modulation, multi-user detection and channel decoding. The example depicted 
here in further detail in Figure 2 is related to equalization. Equalization is re- 
quired because of inter-symbol interference (ISI) resulting from multi-path 
propagation. 

[0024] The equalizer depicted as an example here, comprises a soft interference 
canceller (SC) followed by an MMSE filter (Minimum-Mean- Square Equalizer) 
optimised with the channel decoder feedback, channel response and noise level. 
The SC/MMSE equalizer uses information on transmitted bits. 
[0025] The principle of turbo equalization applied to the usage of SC/MMSE 
algorithm leads to a structure where the likelihood information is exchanged 
between SC/MMSE equalizer and a decoder through the interleaving/de- 
interleaving blocks. The SC/MMSE turbo receiver structure can be applied not 
only the equalization of inter-symbol interference for single-carrier signalling 
but also to various other signalling schemes, such as multiple signal detection in 
multiple-input-multiple-output (MIMO) systems. 

[0026] An example of the SC/MMSE turbo receiver structure is shown in Fig- 
ure 2. The SC/MMSE equalizer 200 is coupled to the channel decoder 204 via a 
de- inter leaver 202. The decoder 204 feeds back a priori information on the 
value of transmitted symbols to the SC/MMSE equalizer via an interleaver 206. 
The method is used in a channel equalizer. The arrows passing between the 
equalizer 200 and decoder 204 through the interleaver 206 and de-interleaver 
202 illustrate the exchange of likelihood information. 

[0027] The decoder 204 can in general be any algorithm capable of providing 
soft symbol information metrics also called symbol decisions, data decisions or 
symbol estimates, for instance. Possible decoder algorithms are, for instance, 
log-domain maximum-a-posteriori (log-MAP), its max-approximation (max- 
log-MAP) and the soft-output Viterbi algorithm (SOVA). 
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[0028] The embodiment of a data processing method in a channel equalizer of a 
receiver utilises the fact that two consecutive covariance matrices called in this 
application a first covariance matrix and a second covariance matrix include 
mostly the same components. The second covariance matrix at the time n+1 is a 
shifted version of the previous (first) covariance matrix at the time n with new 
components only in the last row and column. The approach is especially useful 
in MIMO systems, where multiple desired symbol likelihood information is 
computed with a common interference covariance matrix but different channel 
responses. The equalizer operates with a sliding signal window, and the inter- 
ference cancelled signal covariance matrix can be recursively updated by ex- 
ploiting the knowledge of the memory within the sliding window. The size of a 
sliding step is usually one. 

[0029] Typically, as far as the equalizing process continues, the sliding of the 
signal window continues and the inversion of a new covariance matrix is calcu- 
lated with the aid of the previous inverse covariance matrix utilising common 
covariance components. 

[0030] In cases where the received signal is over-sampled in space or time, the 
covariance matrix has a block structure. The recursion must then be executed 
multiple times per each equalized symbol to match the over-sampling ratio. 
[0031] The embodiment begins in block 300. The embodiment is suitable espe- 
cially for systems where multiple antennas transmit a coded and interleaved sig- 
nal, which is received with one of multiple antennas. In block 302, interference 
in a radio channel is estimated. The interference estimation is known in the art 
and it is therefore not explained here in further detail. One widely used method 
is measuring the impulse response of a received signal. If there are several radio 
channels of interest, the interference is usually estimated in these channels. In a 
multi-user system, the interference can be estimated for every user in a multi- 
user system or a common interference estimate is used. In the example used 
here, the received signal is modelled with the linear model: 
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r = Hb+w, 



(1) 



where 



H is a channel response matrix, 
b is a symbol vector, 



w is a noise vector. 



[0032] On the basis of the interference estimation, a first covariance matrix is 
created. The size of the matrix depends on the amount of receiver sampling and 
the channel response length. The system may be a MIMO (Multiple input, mul- 
tiple output) system, in which case there are several simultaneous radio trans- 
missions to be received through different channels. 

[0033] In block 302, an inverse matrix of the first covariance matrix and its 
Cholesky decomposition are also defined. The inverse matrix of the first covari- 
ance matrix can be calculated using exhaustive calculation used in prior art 
equalizers or another prior art method. 

[0034] The inverse matrix of the first covariance matrix can be partitioned as 
follows: 



Up means a sub-matrix and 

H means a complex-conjugate transpose (Hermitian) matrix. 
[0035] The Cholesky decomposition of the partitioned inverse matrix (2) is de- 
fined as a lower-triangular matrix: 



U(*) = 




P J 



(2) 



where 



Up means a scalar, 
u p means a vector, 

u p H means a complex-conjugate transpose vector, 




(3) 



where 
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cop means a scalar, 
<o p means a vector, 

o H means a zero vector and 

£l p means a lower triangular sub-matrix, 

H means a complex-conjugate transpose (Hermitian) matrix. 
[0036] In block 304, selected covariance components are removed from the 
Cholesky decomposition matrix. The covariance components to be removed are 
typically components related to the observation time of the previous covariance 
matrix. 

[0037] In block 306, the inverse of the sub-matrix E is computed. The sub- 
matrix Z represents the common part of the first covariance matrix and a sec- 
ond covariance matrix, which includes covariance estimates of a second obser- 
vation time, by using the aid of the Cholesky decomposition of the inverse ma- 
trix of the first covariance matrix. 

[0038] Two consecutive covariance matrices called in this application a first 
covariance matrix and a second covariance matrix include mostly the same 
components. The second covariance matrix at the time n+1 is a shifted version 
of the previous (first) covariance matrix at the time n with new components only 
in the last row and column. 

[0039] The sub-matrix £ representing the common part of the two consecutive 
covariance matrices is clarified with the aid of Figure 4. The first covariance 
matrix is marked with number 400 and the second covariance matrix is marked 
with number 402. The sub matrix Z which represents the common part of ma- 
trices 400 and 402 is marked with number 404. 

[0040] The Cholesky factorisation of the inverse matrix Z" 1 can be computed 

as: 

I-' =aa H , (4) 

where 
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SI is a sub-matrix and 

Sl H is a complex-conjugate transpose of the sub-matrix. 
[0041] We find that 

" = (5) 
where 

£l p is a sub-matrix of (3). 

[0042] In block 308, interference from a received signal is estimated at a sec- 
ond observation time and additional covariance components are determined on 
the basis of the estimation. The additional covariance components are the co- 
variance components of the part of the second covariance matrix 402 of Figure 
4 not common with two consecutive matrices, that is to say they do not belong 
to the sub-matrix £ . The location of the additional covariance components is 
marked in Figure 4 with reference number 406. The additional covariance com- 
ponents can be computed for upright part components and then derive the hori- 
zontal part components with the aid of upright components. 
[0043] Hence the additional covariance components can be computed as 
'« 
°J 
where 

c f means covariance component vector, 

a f means covariance component (scalar) (located in the comer of 
the second covariance matrix), 

H(n) means system matrix in the observation time of the second 
covariance matrix, 

diag means a diagonal matrix, 

b(n) is a symbol estimate, 

H means a complex conjugate (Hermitian) matrix, 
<Tq is the noise variance, 
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I is an interference matrix. 



[0044] The computation of new covariance components in (6) is defined for 
BPSK modulation. The extension of the computation to other modulation meth- 
ods can be implemented by known methods. 

[0045] In block 310, Cholesky decomposition of the inverse matrix of the sec- 
ond covariance matrix is created by using unitary rotations. 
[0046] The Cholesky factorisation of the inverse matrix of the second covari- 
ance matrix is defined as a lower-triangular matrix: 



£l f means a lower triangular sub-matrix, 

co" means a vector, 

o" means a zero vector, 
co* f means a scalar, 

H means a complex-conjugate transpose (Hermitian) matrix and 
* means a complex-conjugate. 



[0047] A pre-array is composed with the additional covariance components 
<s f ,a f which are pre-processed with the prior sub-matrix fi. By unitary rota- 
tions, it is possible to iteratively annihilate the last column of the pre-array. If 
the dimension of the matrix W f is d, d-\ rotations are needed to annihilate all 
but the last element of the last column of the pre-array. Thus the Cholesky fac- 
torisation of the inverse matrix of the second covariance matrix can be ex- 
pressed with the pre-array and a concatenated rotation matrix as: 




(7) 



where 




(8) 



where 



SI means a sub-matrix, 
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o H means a zero vector, 

£l H means a complex-conjugate transpose of the sub-matrix, 
u f = (cr f -a^nfi^aj, 

Of means covariance component (scalar), 

0 means series of unitary rotations and 

H means a complex-conjugate transpose (Hermitian) matrix. 
[0048] Algorithms based on unitary rotation are well-suited for pipelined and 
parallel implementations with, for instance, systolic arrays. 
[0049] In block 312, an output value of the channel equalizer is generated by 
utilizing information obtained with the aid of the Cholesky decomposition of the 
inverse matrix of the second covariance matrix. The output value of the channel 
equalizer is typically generated by further utilizing a-priori symbol estimate in- 
formation. The a-priori symbol estimate information may include, for instance, 
symbol decisions (also called, for instance, likelihood information, symbol es- 
timate, data decision). 

[0050] The output of an equalizer can be defined as: 

z k (n) = fi k (n)(a k {n)b k (n) + 7 f ( W ) W/ r (*)), (9) 
where 

<x k (") = 1?(n)tj k (n), (10) 
rj k (n) = W?h k (n), (11) 

A00 = i ° k{n ? _ 2 , (12) 

hi") 

h k (n) is a channel response vector, 
n means an nth symbol, 

H means a complex-conjugate transpose (Hermitian) matrix, 

b k (ri) is a symbol estimate (or a symbol decision, likelihood 

information, symbol estimate or a data decision) based on a channel decoder 
feedback, 
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W f is defined in the form (8), 

?(/*) = r(«)-Hb (13) 
where 

M k (n) is a channel response matrix and 

n means an nth symbol. 
[0051] The computation in (9)-(12) is defined for BPSK modulation. The ex- 
tension of the computation to other modulation methods can be carried out by 
known methods 

[0052] The embodiment ends in block 314. Arrow 316 depicts that the em- 
bodiment may be repeated for calculating the next equalizer output signal value 
usually for the next transmitted symbol. 

[0053] In Figure 5 5 there is depicted a generalised example of an equalizer. It is 
obvious to a person skilled in the art that the equalizer may also include ele- 
ments other than those illustrated in Figure 5. Equalizers are typically part of 
receivers. The receiver may be located in several kinds of communication units, 
such as a base station of a communication system, for instance a base station 
used in a wideband code division multiple access system (WCDMA), such as 
UMTS (universal mobile telecommunications system). 

[0054] Channel matrix H is estimated with a channel estimator and required 
channel response values h k are conveyed from block 502 to block 500 along 
with new (additional) covariance matrix components <* f ,cr f In block 500 the 
Cholesky factorisation of the inverse matrix of the second covariance matrix is 
calculated according to equation (8). The output values of the block 500 are the 
vectors r|D k and W^r where () H means complex-conjugate transpose and is 

r (n) = r(n) - Hb , which means a received signal from where known signal com- 
ponents have been removed. 
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[0055J In multiplier 504, n k («) and w;?(«) are multiplied element-wise. 
Squaring block 506 and adder 510 correspond to the calculation of 7 f («)%(«), 
which gives a k (n) according to equation (10). 

10056] In multiplier 514, a k ( n )is multiplied with b k (n), which is a symbol es- 
timate based on channel decoder 204 feedback depicted in Figure 2. Then the 
result of the multiplication is added to ^W/r(«)in adder 512. 
[0057] In block 518, p k («) is calculated according to equation (12). The output 
of the adder 512 and 0 k («) are multiplied in multiplier 516 and the output of the 
equalizer is then z k («) = p k ( n )(a k {n)b k («) + rj k H («)W/ r («) ) according to equation 
(9). When multiple transmitters are received, h k (n), ^ a k (n) A(#0 and z k (n) 
must be computed separately for each transmitter. 

[0058] Figure 6 shows a generalised example of a receiver. It is obvious to a 
person skilled in the art that the receiver may also include elements other than 
those illustrated in Figure 6. The receiver may be located in several elements of 
a communication system, such as in a base station. In this example, the receiver 
uses the same antenna 600 as a transmitter, and therefore there is also a duplex 
filter 602 to separate transmission and reception. The antenna may be an an- 
tenna array or a single antenna. In the receiver, RF-parts 604 in this case also 
comprise a power amplifier, which amplifies the received signal attenuated on a 
radio path. Typically RF-parts down-convert a signal to an intermediate fre- 
quency and then to a base band frequency or straight to the base band fre- 
quency. The analogue-to-digital converter 606 converts an analogue signal to 
digital form by sampling and quantizing. 

[0059] The equalizer 608 carries out most processing steps of embodiments of 
the data processing method. The structure of the block is described in Figure 5 
in general. Arrow 612 depicts that symbol estimates from a channel decoder are 
fed to equalizer 608. In this embodiment, the channel decoder is a part of DSP 
block 610. Typical functions of a DSP-block are, for example, de-scrambling, 
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decoding, de-interleaving, etc. The basics of the digital signal processing are 
known in the art. The equalizer may also include means for channel estimation. 
The interleaver and de-interleaver are not depicted in Figure 6 for the sake of 
clarity. 

[0060] The embodiments of the data processing method are typically imple- 
mented as a processor and software, but different hardware implementations are 
also feasible, e.g. a circuit built of separate logics components or one or more 
client-specific integrated circuits (Application-Specific Integrated Circuit, 
ASIC). A hybrid of these implementations is also feasible. Additionally, unitary 
rotations can be implemented with systolic arrays. 

[0061] Even though the invention is described above with reference to an ex- 
ample according to the accompanying drawings, it is clear that the invention is 
not restricted thereto but it can be modified in several ways within the scope of 
the appended claims. 



